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Abstract 

We discuss D-dimensional Euclidean scalar field interacting with a 
scale invariant quantized metric. We assume that the metric depends on 
d-dimensional coordinates where d < D. We show that the interacting 
quantum fields have more regular short distance behaviour than the free 
fields. A model of a Gaussian metric is discussed in detail. In particular, 
in the "l? 4 theory in four dimensions we obtain explicit lower and upper 
bounds for each term of the perturbation series. It comes out that there is 
no charge renormalization in <1> 4 model in four dimensions. We show that 
in a particular range of the scale dimension there are models in D = 4 
without any divergencies. 

1 Introduction 

We discuss D-dimensional Euclidean scalar field interacting with a quantized 
scale invariant metric. The metric depends on d-dimensional coordinates. The 
simplest case, which arises from a Gaussian metric, will be discussed in detail. 
If d = 2 then we can give at least two non- Gaussian examples of such a scale 
invariant metric. In one example the metric is described as a two-dimensional 
SL(D, R)-va\ued field. As the second example we may consider the Polyakov 
model [Jy in a non-critical dimension which in the conformal gauge is reduced 
to the two-dimensional Liouvillc model. 

In the first model we treat a metric tensor on the Ricmannian manifold 
(Euclidean formulation) 

(G) AB = g AB (1) 

as a two-dimensional field G with values in a set of real symmetric positive 
definite D x D matrices G. We choose the metric in a block diagonal form 
g ab = S AB if A,B > D -2 and for A, B < D - 2 the tensor /"(xy) is a 
(D — 2) x (D — 2) matrix depending on xj? 6 R 2 . The manifold of positive 
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definite matrices is homeomorphic to B x SL(D — 2, R)/0{D — 2). We choose 
a conformal invariant action for G 

W{G) =TrJ d* F G- l dGG- l dG + WZW (2) 

where d = d\ — idi is the holomorphic derivative and WZW denotes the Wess- 
Zumino-Witten term Q . The conformal invariance of this model has been shown 
inrefs. [1|§. 

The first model is not invariant under the group of diffeomorphisms of the 
metric. In the second conformal invariant model of gravity we consider the string 
coordinates interacting with two-dimensional gravity in a way invariant 
under general coordinate transformations 

W = J dx F Jgg ab d a X»d b X» (3) 

The classical model does not depend on the metric but the quantum one depends 
on the conformal Weyl factor. We choose g a b — Sab exp \ ■ A functional integral 
over X leads to the effective action 

w eff (g) = J dx F (OaXdaX + <* exp x) (4) 

This Liouville model is scale invariant || . 

In sees. 2-3 we discuss general scale invariant models. In sec. 4 we restrict 
ourselves to the Gaussian metric where we can obtain detailed estimates on the 
perturbation series of the (<&<!>* ) 2 interaction ( we discuss the complex scalar 
field instead of the real one just for simplicity of the Gaussian combinatorics). 
It is shown that quantum fields interacting with a singular random metric are 
more regular than the free fields (a conjecture reviewed in [||; see also 0) 



2 The scalar propagator 

We consider a complex scalar matter field $ in D dimensions interacting with 
gravitons depending only on a c?-dimensional submanifold. We split the coor- 
dinates as x — (xc,Xi?) with xp £ R d . Without a self-interaction the 
correlation function is equal to an average 

ft J Vgex V ^W{G)^A- 1 {x iy ) (5) 

over the gravitational field g of the Green's function of the operator 

D-d D 

- A =2 ^ <T(x F HA + - d l 

fi=l,v=l k=D-d+l 
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We repeat some steps of ref.[[lO| (our case here is simpler and more explicit). 
We represent the Green's function by means of the proper time method 



A 1 (x,y) = / dT(exp(-TA))(x,y) (7) 
Jo 

For a calculation of (exp (— tA)) (x, y) we apply the functional integral 

K T (x,y) = (ex P (-rA))(x,y) = /©xexp(-i / - \ f g^)******) 

5 (x (0) - x) 6 {x (r) - y) 

In the functional integral (8) we make a change of variables (x — > b) determined 
by Stratonovitch stochastic differential equations Jll| 

dx u (s)=e%(x(s))db A (s) (9) 

where for Cl = 1, 2, D — d 

and = 5% if > D - d. 

As a result of the transformation x — > 6 the functional integral becomes 
Gaussian with the covariance 

E[b a {t)b c (s)} = S ac min(s,t) (10) 

In contradistinction to [^0) eq.(9) can be solved explicitly. The solution q T of 
eq.(9) consists of two vectors (qcqf) where 

q F (r,x F ) = x F + b F (r) (11) 

and qc has the components (for fi = 1, D — d) 

q»(r, x) = + f e£ (q F (s, x F )) d6 a ( S ) (12) 
Jo 



The kernel is 



^ T (x,y)= J B[%-« T (x))] = 

= £[<5(y F - x F - b F (r)) ^ J (y„ - q^ (r, x))] 



(13) 



Using eq.(12) and the Fourier representation of the 5-function we write eq.(13) 
in the form 

K T (x,y) = (2ir)- D f dp G dp F 

E[exp (ip F (y F - x F ) + ip G (y G - x G ) - ip F b F (r) - i J (q (s, x F )) db a (s))] 

(14) 
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3 The scale invariant model 

In general, we cannot calculate the average over the metric explicitly. How- 
ever, the scale invariance of the metric is sufficient for a derivation of the short 
distance behaviour of the scalar propagator. 

Let us note that y/rb(s/T) ~ b(s) where b denotes an equivalent Brownian 
motion (the equivalence means that both random variables have the same cor- 
relation functions). Then, using the scale invariance of e with the index 7 we 
can write 

e(Vrx F ) ~ r~i e(x F ) (15) 

Hence, in eq.(12) 

g"(T,x) = x» + t5-3 J g£ (j-ixF + bp (*)) db a (s 

The expectation value over e is 
(K T (x,y)) 



(16) 



-^(1-7)- 



■ (E [d ((y F - x F ) T -i - b F (1)) S (r-i+i (y - x) - 77)] ) 

(17) 

where 



if = J ~e» (r" W + b F (*)) db a (s) 



Let P(u, v) be the joint distribution of (r/,b F (l)) ( P does not depend on 
x F because of the translational invariance). Then, the propagator of the $ field 
is 



h{A-\x,y)) = 

hj °° dTT-V-^-tP ((x G - y G ) r(-^)/ 2 , (x F - y F ) r^) 



(18) 



Eq.(18) in momentum space has the representation 

f°° ~ 1-1 

hiA-^kckp)) =h dTP(T—k G ,V^k F ) (19) 
Jo 

where P denotes the Fourier transform of P. Using eq.(14) we may write 

= h dr(E[expi (yrkpbp (1) + k G ?7 G )]} 

The dispersion relation (relating the frequency to the wave number) is deter- 
mined by (after an analytic continuation k — > ifco) 

((^l- 1 (k G ,k F )»- 1 =0 
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It can be concluded from eq.(18) that in general the dispersion relation will 
be different from the standard one (resulting from a wave equation) ko ~ |k|. 
In particular, we can see that if |k F | > |k G | then (^(kckf)) - IM" 2 
whereas if |k G | S> |kp| then (.4 _1 (k G , k F )) ~ |k G |~~. In the configuration 
space, the propagator tends to infinity if both |x F — yp\ and |x G — yd tend 
to zero. However, the singularity depends in a rather complicated way on the 
approach to zero. It becomes simple if either |x F — y F | — or |x G — yd =0 . 
So, if |xf? — yp\ = then we make a change of the time variable 

r = i|x G -y G |T^ (20) 

Using eq.(18) we obtain the factor depending on |x G — y G | in front of the integral 
and a bounded function A of coordinates ,i.e., 

(A- 1 (x,y))=A\x G -y G \- D + 2 

If |x G — yd =0 then we change the time variable 

t = t\x F - y F \ 2 

As a result 

{A-\x,y)) = A|x F -y F r (D - 2)(1 - 7) (21) 

with a certain bounded function A. We can see that in the x G coordinate the 
singularity remains unchanged but the propagator is more regular in the x F 
coordinate. 

It is not possible to calculate the probability distribution P exactly. Choos- 
ing as a first approximation r/ ~ b G (l) we obtain 

P(u,v) = (2 7 r)-*cxp(-^-^) 

In this approximation 

ft r°° 1 1 

h(A- 1 (k G ,k F )) = ^J o drexp(--r 1 ^|k G | 2 --r|k F | 2 ) (22) 

In the Higgs model we need the mass term j : $$* :. The perturbation series 
for the mass will be finite if the integral 

/ dxpclxoKA-^x^y))] 2 
Jv 

is finite for any bounded region V. This integral is convergent if 

/ dkGdkFKA^ik)}] 2 < oo (23) 

J\k\>e 
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for any e > 0. It can be checked that if D = 4 then for any 7 > the integral 
(23) is convergent. 

In fact, in D dimensions the convergence of the integral (23) (with the proper 
time representation (18)) follows from the convergence of the integral 



J Jo 



(24) 



It is finite if d + (1 — j)(D — d) < 4. Together with the inequality 7 < i wc 
obtain the inequality D < 8 — d. 

The integral (24) comes directly from the approximation (22) . However, we 
can prove this result in general under some mild regularity assumptions on the 
Fourier transform (19) of P defined in eq.(18). So, 

J \{A- 1 (k))\ 2 dk = J drdr'P(r^k G ,\/7k i ,)P(T' i ^ 1 k G ,v / 7 7 k F ) ( ik G dk F 

We introduce the spherical coordinates on the plane r = rcos^,r' = rsin#. 
Then, just by scaling of r we derive the result (23) under the assumption that 
the integral over 6 is finite. 

4 The ($<r) 2 model in a random Gaussian met- 
ric 

After the general scale invariant models of gravity in the previous sections we 
consider now a Gaussian model. For the Gaussian model we can prove explicit 
upper and lower bounds on the correlation functions. We consider a complex 
scalar matter field $ in D dimensions . Using eq.(12) and the Fourier represen- 
tation of the ^-function we write eq.(14) in the form 

K T (x, y) = (2TT)- D+d f dp G exp (ip G (y G - x G )) 

E[S (y F - x F - b F (r)) exp (-» J Pfi e% (q (s, x F )) db a (s))} 

The random variables b F and b a are independent. Hence, using the formula 

£[expz j f a (q F )db a ] =£[exp(-i J fafads)} 
we can rewrite eq.(14) solely in terms of the metric tensor 



K r (x, y) = (2vr) D+d J dp G exp (ip G (y G - x G )) 

E[S (y F - x F - b F (t)) exp (-| / p^g^ (q (s, x F )) ds)\ 



(25) 



We assume that the metric g^ v is Gaussian with the short distance correlations 
(g^(x F ) 5 ^(y F )) - -DW>(x F - y F ) = -C™|y F - x F |- 4 ^ (26) 
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where C is a scale invariant tensor and D must be positive definite if the momen- 
tum integrals in the final formula are to exist. Such a requirement contradicts 
the positive definiteness of the action for the gravitational field. However, in the 
Einstein gravity the conformal modes give a negative contribution to the action. 
The model of a conformally flat metric with C pv '' pp = S plJ S crp would be satisfac- 
tory for our purposes ( a proper contour rotation in the complex space of metrics 
is needed in order to perform the functional integral; such an interpretation of 
the functional integral over the conformal modes has been considered also in 
quantum gravity p^||l3| ). Then, the integral (17) over g can be calculated 

(K T (x, y)) = (2tt)~ d J dp G exp (ip G (y G - x G )) 

E[6 (y F - x F - b F (r)) exp (-1 f Q T V»V*VvV P D^ ap (b F (s) - b F («')) dsds')} 

(27) 

The propagator of eq.(3) has the form 

(A- 1 (x,y)) = / drT-i^ D - d ^-^/ 2 F(T-i(y F - XF ),r-i+i(y G -^ G )) 
Jo 

(28) 

There is a restriction on the allowed singularity of the two-point function of the 
metric field. So, if the random variable in the exponential (25) is to be well- 
defined (without any renormalization) then the expectation value of its square 
should be finite ^ 

(E[( T g^(q( S) x F ))ds) ]) < oo 



This expectation value leads to the integral (if the singularity of D is |u F | 47 ) 

J du F J dsj^ ds'(s - s')~= exp ( - iu F /(s-s'))|uFr 47 

The integral is finite if 7 < \. 

We calculate the expectation value in the gravitational field (26) of a product 
of any number of propagators. First, consider the two-point function of the Wick 
square which we need for the mass term in the Higgs model 

(: ^(x) : (y)) = ((A- 1 (x,y)) 2 ) 

= (2tt) 2D+2d J dTidr 2 J dp G dp' G exp (ip G (y G - x G ) + ip' G (y G - x G )) 
E[S (y F - x F - b F (n)) S (y F -x F -b' F (r 2 )) 

exp ( - \ Q / ri Vp.VaVvVpD^ ap (b F (s) - b F (s')) dsds' 

-I Jo" Jo 2 P>'M P D^ P (b' F (s) - W F (s>)) dsds' 

~\ Jo 1 Jo* P^P' P P',D^ P (b F (s) - b' F (sO) dsds')} 

(29) 

We introduce the spherical coordinates on the (r l7 r 2 )-plane T\ — rcosQ and 
T2 = r sin#. Next, we rescale the momenta k F = p F \Zr, k.' F — p' F \fr , k G = 
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Par 2 % and k G = p' G r 2 % . Then, we can see that 
(: $$* : (x) : <M>* : (y)) 

= / dOdrrr-^-^-VFie, r~b (x F - y F ), r^+i (x G - y G )) 



(30) 



It follows just by scaling of coordinates (the r- integral scales as twice the r- 
integral in eq.(18)) that for short distances 

(: $$* : (x) : $$* : (y)) ~ ((^^ (x,y))) 2 (31) 

We need to prove that correlation functions (28) are finite and non-zero. We 
show first that the bilinear form (fj, [A~ 1 )fi) is finite and non-zero on a dense 
set of functions /. For this purpose we choose 

/k(x G ) = (27ra) _£ ^ exp(-^x G + ikx G ) 

Then, 

{hM'^hf) = (2K)- D+d J™ drr~i jdp G 

E[S fr-i(y F -x) F -b F (l)) 

exp(-^(p G -k) 2 -i(p G -k') 2 (32) 

-\^ 21 lo P»P«P»PpD m (b F (s) - b F ( s ')y s ds> 

Both sides depend on x F and y F because we integrated out only x G and y G . 
In our estimates we apply Jensen inequalities in the form (for real functions A 
and /) 

E[exp A] > exp E[A] (33) 

and 



£[exp ( 



1 , /•! 



dsds'f(s lS ')^}< dsds'E[exp{-f(s,s'))} (34) 



An upper bound can be obtained by means of the Jensen inequality (34) ex- 
pressed in the form 

< 2 J Q dr J Q ds J ds' J du.1du.2dpc 

r-i exp (-i (p G - k) 2 - i (p G - k') 2 ) p(s', ui)p(* - s', u 2 - ui) ( 35 ) 
P ( 1 - s, (y - x) - u 2 ) exp ^- z ^ 1 PpPaPuPpD^ ap (m - u 2 )J 

where p(s, u) = (2tts)~^ exp(— u 2 /2s). We can convince ourselves by means of 
explicit calculations (using a proper change of variables) that the integral on 
the r.h.s. of eq.(35) is finite. For the lower bound it will be useful to introduce 
the Brownian bridge starting from x and ending in x + u [ Oj defined on the 
time interval [0, 1] 

a(x, u, s) = x + us + c(s) (36) 



8 



where c is the Gaussian process starting from and ending in with the corre- 
lation function 

E[c j (a')ck(a)] = 5 jk a'(l-8) 

for s' < s. Then, the 6 function in eq.(27) defines the Brownian bridge and the 
Jensen inequality (33) takes the form 

(/k.M" 1 )/*) > (2n)- D + d J co drr-ijdp G 

eX P ( ~ h (PG ~ ~ h (PG - k ') 2 - jT 2 -^ jl PuPaPvPp 

E[D^ a P (a (o, T-iy F - t-5x f , s) - a (o, r"iy F - r-5x F ,s'j )]dsds') 

(37) 

where the expectation value in the exponential on the r.h.s. of eq.(37) is equal 
to 

J du J ds ds'uj(s, s')~i exp ^ — iu 2 /w(s, s')^ 
|u - (y F - x F ) + T-h' (y F - x F ) |~ 47 



(38) 



where w(s, s') = (s - s'){l - s + s'). It is finite if 7 < \ (the form (26) of the 
graviton two-point function is assumed). 

We can make the same estimates for the Wick square 

(: $$* : (/ k ) : : 

= (27r)- 2 ^+ 2d / dndr 2 J dp G dp' G exp (- £ (p G + p' G - k) 2 - £ (p G + p' G - k') 2 ) 

^[^ (YF - x F - b F (n)) 5 (y F - x F - b' F (r 2 )) 

exp ( - i J T1 J T1 p^PuPpD^P (b F (s) - b F (*')) dsds' 

4 Jo" Jo" (b' F (*) - b' F (*')) 

"§ /eT Jo" ?wX^^ P (b F (*) - b' F (*')) <w)] 

(39) 

It is clear that we can obtain lower and upper bounds applying the Jensen 
inequalities (33) and (34). So, for the upper bound 

(: $$* : (/ k ) : <S»S>* : (.f k ')> 

< 2{2ir)- 2D+2d J d Tl dT 2 Jl ds J Q S ds' / dp G dp' G 
exp (-£ ( PG + p' G k) 2 - £ (p G + p' G - k') 2 ) 
E[S (y F - x F - Vnbf (1)) S (y F - x F - V^ b ' F C 1 )) (40) 
exp ( - ^WW^^ (^b F (s) - Vnb F («')) 
-\r^p' a p' v p' p D^P (V^b F (a) - V^ b f (*')) 
-Ir^p^p'^D^P (Vnb F (a) - ^b' F (s')) 



where the r.h.s. of cq.(40) can be expressed by the transition function for the 
Brownian motion as in eq.(35). 
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For the lower bound we obtain 

(: *$* : (/ k ) : : (/ k ,)) > (2n)- 2D + 2d J d Tl dr 2 J dp G dp' G 
r^Kp cxp (-£ (p G + p' G - k) 2 - X ( PG + p G _ k ') 2 ) 
cxp ( - E[\ ft PvPaPvPpDWP (a (a) - a {s')) dsds' (41) 
-I Jo JoP»lP' p D^ (a' (,s) - a' (s')) dsds' 
\ So Jo' *WX^ ;CT " (a (*) - a' (*')) dsds']) 

where we denoted 

a(s) = x F + (y F - x F )s + v^Tlc(s) 

and 

a'(s) = x F + (y F - x F )s + v^c'(s) 

The r.h.s. of eq.(41) can be calculated explicitly using the correlation function 
for the Brownian bridge. 

We compute now higher order correlation functions 

(<f>(x)$(x')$*(y)$*(y')) 

= (A- 1 (x,y)A- 1 (x',y')) + (x^x') 

= (2n)- 2D+2d J d Tl dT 2 J dp G dp' G exp (ip G (y G - x G ) + ip' G (y' G - x' G )) 
E[S (y F - x F - b F (n)) 5 (y F - x' F - b' F (r 2 )) 

exp ( - \ / T1 J T1 Plx p Mp D^P (b F (s) - b F (a')) dsds' 

-3 JcT Jo" P '^>y p D^" (b' F (s) b' F (s>)) dsds' 

"2 Jo" S? P^P' P P'.D^ P (x F - x F + b F (s) - b' F (s')) dsds') ] + (x -> x') 

(42) 

where (x — > x') means the same expression in which x is exchanged with x' . 
The four-linear form (42) calculated on the basis / reads 

($(/ kl )$(/ k3 )$*(/k 2 )**(/k 4 )) 

= (27r)- 2B + 2d / dndT 2 J dp G dp' G E[6 (y F - x F - b F ( Tl ))6 (y' F -* F -b' F (r 2 )) 
exp (-£(p G - kO 2 - ^(p G - k 2 ) 2 - £(p' G - k 3 ) 2 - £(p' G - k 4 ) 2 ) 
exp ( - \ J T1 J T1 PflPMp D^P (b F (s) - b F (s')) dsds' 
-I Jo" Jo" WM^™ { y p {s) Wp («')) dsds' 

-I Jo" Jo' 2 P^P'pP'.D^ (x F - x F + b F (s) - b' F (s')) dsds')] + (i, 2 -> 3, 4) 

(43) 

where the last term means the same expression with exchanged wave numbers. 
We introduce the spherical coordinates on the (n, T 2 )-plane n = rcosO and 
r 2 = rsin0. Let us rescale the momenta k = P\/r, k' = p'y 7 ^ , k = pr?~2 
and k' = p G r 2 ~ 2 . Then, we can see that the four-point function (42) takes the 
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form 

($(x)$(x')$*(y)$*(y')) 

= f d6drrr- d -^-^ D -^F 4 (e,r-i(x F - y F ),r-i{x' F - y' F ), 
r-5(x^-x F ),r-5(x' F -y F ),r-2(y^ -x), 

r-5 + i( XG -y G ),r-3 + i( x ' G -y^),r-3 + i( x ' G - y G ), r -3 + i ( XG - y' G )) 

(44) 

We can conclude from eq.(44) just by scaling that the singularity of the four- 
point function is a product of singularities. 

We must define now the J($$*) 2 interaction . First, let us calculate the 
two-point function of the interaction Lagrangian 

(: :(x):($$*l:(y)) = ((A- 1 (x,y)) 4 ) 

= (27r)- 2D + 2d j E[U a a Zt dr a d V a exp (i V a (y G - x G )) 6 (y F - x F - b a F (r„)) 
exp ( - | E„,.- Jo" So" Vlv%viviv^° p (b% (s) b F ' (*')) dsds')] 

(45) 

For the lower bound it is sufficient if we let x G = y G . Then, 
(: (<M>*) 2 : (x) : (<M>*) 2 : (y)) 

> (27r)- 2 ^+ 2d / n::t dr a d P <v- * (46) 

exp ( - |SE„,„/ Jo IoP»>p' D ^' ap ( a F («) - 4 (*')) <*"**')] 

where the expectation value in the exponential can be calculated and the integral 
on the r.h.s. is finite. For the upper bound of the interaction J v dx($$*) 2 (x) 

we may take the test functions of eq.(32) with 2ira = V^^ 1 and k = in order 
to approximate the finite volume integral. With these test-functions we obtain 
the upper bounds in the same way as we did it in eqs.(35) and(40) from the 
Jensen inequality (34). 

Through an introduction of spherical coordinates in the (n, T4) space we 
can show that for short distances 

(: ($$*) 2 : (*) : (d>$*) 2 : (y)) ~ ((A^ (x,y))) 4 (47) 

(because the r-integral scales as four-times the T-integral in eq.(28)). 
Let us calculate now the vacuum diagram 

h = f dx f dy(: (<i><I>*) 2 ■ (x) : (4><I>*) 2 : (y)) (48) 
Jv Jv 

corresponding to the second order perturbation expansion in the coupling con- 
stant. In the momentum space the convergence of this diagram follows from the 
convergence of the integral (for large momenta) 

J dk 1 dk2dk 3 {A- 1 (k 1 )A- 1 (k2)A- 1 (k 3 )A- 1 (k 1 + k 2 + fc 3 )> (49) 
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On the basis of eq. (47) the integral (49) can be approximated by 

J dhdhdksiA-^k^iA-'ik^iA-^hViA-'ih + ki + h)) (50) 

In the proper time representation the convergence of the integral (50) depends 
on the convergence of the integral (for small r) 

d — i 

J dT 1 dT 2 dT 3 dT 4 (T 1 + T 4 ) 2 (-TfrJ + (TlT 3 + TlT 4 + T3T4) (T1T2 + T1T4 + T 2 T 4 )) 2 

(t 1 1 - 7 +t 4 1 -' 1, )-*+t 

( - rr^rr 2 ^ 

+ (r 1 1 -V 3 1 ^ + tv 1 "^ + r 3 1 ^r 4 1 ^)(r 1 1 -'r 2 1 ^ + r^r^ + r^r^)) l_ " 

(51) 

We obtained the formula (51) from the representation (19) by scaling of mo- 
menta as in the argument at the end of sec. 3 (we could obtain lower and upper 
bounds on the expectation value (48) using the lower and upper bounds on the 
Wick powers as in eqs.(40) and (46); however the explicit integrals are harder 
for analysis than eq.(51)). In order to investigate the convergence of the integral 
(51) it is useful to introduce the spherical coordinates. Then, the condition for 
the convergence of the radial part of the integral (51) reads 

|d+|(D-d)(l-7)<4 (52) 

For D = 4 this condition has a solution only if d = 1 and 7 > | (together with 
7 < i). If the condition (52) is satisfied then we obtain quantum field theory 
without any divergencies. However, with any 7 > the ultraviolet behaviour 
is better than the one with 7 = 0. As can be deduced from the dimensional 
regularization, if a diagram is only logarithmically divergent (as for the charge 
renormalization in D = 4 model), then after coupling to quantum gravity it 
becomes convergent. In D = 4 with 7 > (but without the inequality (52)) 
the vacuum diagrams are divergent but the ones corresponding to the charge 
renormalization are convergent. We calculate, e.g., the four-point function in 
perturbation expansion at the second order in the coupling constant 

J d i z l d i z 2 ^{x l )^{x 2 )^{y l )^{y 2 ) : ($$*) 2 : ( Zl ) : (<&<T) 2 : (z 2 )> (53) 

It is clear from our estimates that the four-point function (53) will be non-trivial 
and finite with any 7 > 0. We can calculate n-point functions. It follows that 
with the lower and upper bounds established in this section we can obtain finite 
lower and upper bounds on each term (eventually after a mass renormalization) 
of the perturbation series of the (<&<!>*) 2 model. 
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The Gaussian metric with the two-point function (26) may be unphysical. 
We applied this model in order to obtain explicit estimates on correlation func- 
tions in the perturbation series. As shown in sec. 3 the regularizing property 
of the singular random metric is universal. We expect that with some harder 
work the estimates on the perturbation expansion in ($$* ) 2 , with the two- 
dimensional gravity of sec.l, are possible as well. 
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